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Abstract 

Background: Ihe etiology of birth defects lias been widely studied but is not yet fully clarified, previously published data 
had suggested that maternal age or parity maybe involved, but without consistent conclusions. 

Methods: A population-based, case-control study was nested in a cohort of perinatal infants born from 2010 to 2012 in 
Baoan District, Shenzhen. Four categories of isolated birth defects were defined as cases: congenital heart defects (CHD, 
n = 693), Polydactyly (n = 352), cleft lip with or without palate (CL/P, n = 1 59) and equinovarus (n = 1 1 9). Controls were non- 
malformed infants (n = 1 1 ,307) randomly selected from the same area and period. Odds ratios (ORs) and the 95% confidence 
intervals (CIs) were computed by multivariable unconditional logistic regression analysis. 

Resu/ts: Young maternal age {<25 years old) was associated with a reduced risk of CHD (adjusted OR = 0.73, 95% CI 0.59- 
0.90), while with an elevated risk of Polydactyly (adjusted OR =1.42, 95% CI 1.09-1.84). Increased risk of CL/P-affected 
pregnancy was observed in mothers older than 35 years old (adjusted OR = 2.1 2, 95% CI 1.26-3.57). Compared to primipara, 
those having their second, and third or more delivery were less likely to have infants with equinovarus, with significant 
adjusted ORs of 0.59 (0.40-0.89) and 0.42 (0.19-0.93), respectively. 

Conclusion: IVlaternal age was significantly associated with CHD, Polydactyly and CL/P relevant pregnancy. Mothers with 
higher parity might have lower risk of equinovarus occurrence in offsprings. 
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introduction 

The etiology of non-chromosomal birth defects has been widely 
discussed but is not yet fully clarified, though numerous studies 
have confirmed that environmental and behaviors components are 
involved [1-6]. Apart from these influences, previous studies have 
implicated that maternal inherent characteristics, for instance, age 
[7-10] or parity [9-13] might be also responsible for certain 
categories of congenital defects. However, no consensus was 
reached, and results were usually inconsistent due to limited 
sample size and methodological difficulties. 

Baoan district is accounted for the second largest area among 
the six functional districts in Shenzhen City, and accounted for the 
largest proportion (38.8%) of Shenzhen's total population 
(10.36 million) at the end of 2010 [14]. The average birth rate 
(15.88%o) and incidence rate of birth defects (14.14%o) were also 
reported higher than those of other districts [14]. 

As is generally accepted, population characteristics differ 
considerably across geographical location, therefore public health 
decision making should be based on the specific regional 
researches. However, etiologic studies regarding congenital defects 



in Chinese southern area have rarely been reported in recent 
years. 

In the light of the above, we conducted a case-control study to 
investigate the relationship between maternal age, parity and 
selected isolated birth defects through a large cohort from Baoan 
district, Shenzhen, China. 

Metliods 

Study Design 

This study is a multi-center, population-based case-control 
study within a cohort of perinatal infants born in Baoan district 
from January 1", 2010 to December 3V\ 2012. Data utilized in 
this study were derived from Shenzhen Maternal and Child 
Health Management System (available on http://218.18.233. 
230:8888/), a continuously updated health care electronic 
database with an approximate 99.5% geographic coverage of 
the total area of Shenzhen. Anonymity and confidentiality were 
guaranteed, and written consents were given from the expectant 
parents, or their next of kin on the behalf of the newborn 
participants, for their information to be stored in the System at the 
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time of the expectant mothers were awaiting delivery in the 
hospital. The study protocol was approved by the institutional 
review board of Baoan Maternal and Child Health Hospital. 

Case and Control Ascertainments 

Perinatal infants (from 28 weeks' gestation to one week after 
birth) with congenital malformations were identified by Birth 
Defects Monitoring System (BDMS), a branch of Shenzhen 
Maternal and Child Health Management System. 

In the present study, cases were defined as perinatal infants 
confirmed with isolated congenital heart defects (CHD), Polydac- 
tyly, cleft lip with or without palate (CL/P) or equinovarus, which 
represented the top four phenotypes of birth defects in Baoan 
district. Those with multiple congenital defects, syndromic diseases 
or chromosomal aberrations were excluded as high genetic 
components were concerned [5]. 

The control pool was composed of all non-malformed perinatal 
infants bom in the area from which the cases were recruited 
during the same period (2010-2012). Out of 203,154 records that 
met the primary criteria, 189,378 (93.2%) with complete 
information. The Power and Sample Size Calculations program 
(PS, version 3.0.43) was used to assess the minimum sample size of 
controls to guarantee a power of at least 79.6% (alpha 0.05, two- 
sided) to detect an OR of 2.0 or greater. Eventually a number of 
11,305 records were randomly selected and constituted a share 
control group. 

Data Collection 

Data on maternal and child records were prospectively entered 
into the Shenzhen Maternal and Child Health Management 
System. Once the cases and controls were ascertained, their basic 
characteristics were retrospectively reviewed and the selected 
covariates were extracted independently by two authors (Luo and 
Gao): general information of mothers (nationality, household 
registry, etc) and prenatal infants (gender, gestational age, birth 
weight, etc), the key exposure variables (maternal age and parity) 
as well as other covariates that were chosen to included in the 
analysis. Once the data extraction was completed, the outliers 
were rechecked by looking through the original medical records. 

With rapid industrialization and urbanization, millions of 
Chinese peasants left their hometowns and headed to Shenzhen 
for job, making this city a famous modern industrial zone and 
worker-gathered area. For better management, Shenzhen govern- 
ment has classified the inhabitants' household registry into three 
types: 1) local resident for those who achieved permanent 
residence permit, 2) temporary population for those who lived in 
Shenzhen for more than one year without permanent residency 
and 3) migrant population for those who lived in Shenzhen for less 
than one year without permanent residency. The frequency of 
antenatal checkup was defined as the total numbers that had been 
taken through the whole gestation and treated as a categorical 
variable with four levels: 0, 1-5, 6-10 and >1(). We al.so further 
grouped the checkup into "regular" and "irregular". Regular 
checkup was defined as the first examination be taken from the 
12**" weeks of gestation and more than five times be counted 
through the whole gestation. 

In this study, maternal age was categorized into four groups: 
<25, 25-29, 30-34 and S35 years old. The second group was 
chosen as a reference because generally women of 25-29 years old 
has the lowest overall prevalence and represents the largest 
proportion. Parity was defined as the numbers of live birth a 
female has borne prior to the index pregnancy [15], and 
categorized into three levels: 0 (primipara, reference), 1 and &2. 



Statistical Analysis 

For the enumeration data, independent sample test or 
Fisher's exact test were performed to compare the status of case 
and control groups. Multivariable unconditional logistic regression 
analysis was used to compute the odds ratios (ORs) and the 
corresponding 95% confidence intervals (CIs) for the effects of 
maternal age and parity on the selected malformations. Two 
regression models were used: one crude model without any 
adjustments and another adjusted model that included adjust- 
ments for covariates that were deemed to be plausible confounders 
(covariates that were differentially distributed between cases and 
controls). Specially, liner regression analysis using a stepwise 
procedure was carried out to test co-hnearity problems (eg. 
gestational age and birth weight) prior to multivariable logistic 
analysis. Variance inflation factor (VIF) SlO indicated strong co- 
linearity and then only one significant covariatc- was sck-cted into 
the adjusted model. All statistical analyses were conducted by the 
SPSS, version 13.0 (Chicago, IL, USA). P value <0.05 was 
considered statistical significance. 

Results 

Cases and Controls Enrollments 

A number of 693 CHD, 352 Polydactyly, 159 CL/P and 119 
equinovarus cases and 11,307 unaffected healthy controls were 
enrolled in the final analysis. The selected information of study 
subjects were summarized in Table 1. 

In many instances, there were statistically significant differences 
between case and control groups. A male predominance was noted 
among the Polydactyly (64.8%) and CL/P (64.2%) cases as 
compared with controls (54.1%, P<0.05). Infants with CHD and 
CL/P in general had a lower birth weight and were more often 
prematurely bom than those without malformation (P<0.01). The 
frequencies of still birth were also statistic;Jly higher in the CHD 
(8.2%), CL/P (15.1 %) and equinovarus (2.5 %) groups than those in 
the control group (0.3%). As the liner regression results had obviated 
the possibility of strong co-linearity (VIF= 1.439) between gesta- 
tional age and birth weight, the above and other significantly 
incomparable covariates were considered to be biologically 
potential confounders and were entered into the adjusted model. 

As listed in Table 2, more CL/P or erjuinovarus case-mothers 
fell into the "migrant population" group in comparison with 
control-mothers (76.7% vs. 61.7%, 74.8% vs. 61.7%, respectively). 
Significant lower frequency of antenatal checkup was also seen in 
CL/P-afiFected mothers when using 5 times as cut-off (85.5% vs. 
68.6%). On the other hand, matemal pregnancy complications 
such as gestational diabetes mellitus or hypertension were sparse, 
thus most of the associations were detected as non-significant. 

Co-linearity test for maternal age and parity 

As there exists a possibility that multi-parity is simply acting as a 
marker for maternal age, we also undertook liner regression 
analysis to test co-linearity problems, which were not demonstrat- 
ed by the test results (\TF= 1.184), indicating that the impacts of 
parity were independent by maternal age. 

Maternal age and risk of selected isolated birth defects 

Table 3-4 showed the risk estimates of maternal age, parity in 
relation to the risk of delivering an infant with certain isolated 
birth defects. A protective effect on the occurrence of CHD was 
seen in younger mothers (adjusted OR = 0.73, 95% CI 0.59-0.90). 
On the contrar)', obser\'ations of the Polydactyly group indicated 
that mothers younger than 25 years old were more likely to have 
infants with Polydactyly, albeit with modest association (adjusted 
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Table 1. Basic characteristic of cases with isolated birth defects and unaffected controls. 





Characteristics 


CHD " 


Polydactyly 


CL/P " 


Equinovarus 


Control 


(n = 693,%) 


(n = 352,%) 


(n = 1 59, %) 


(n = 1 1 9, %) 


(n = 1 1 307, %) 


Gender ** 


Male 


396 (57.1) 


228 (64.8) 


102 (64.2) 


59 (49.6) 


6117 (54.1) 


Female 


297 (42.9) 


124 (35.2) 


56 (35.2) 


60 (50.4) 


5190 (45.9) 


Gestational age (week) 


<32 


46 (6.6) 


3 (0.9) 


11 (6.9) 


2 (1.7) 


72 (0.6) 


32-37 


185 (26.7) 


52 (14.8) 


38 (23.9) 


21 (17.6) 


1349 (11.9) 


2:38 


462 (66.7) 


297 (84.4) 


110 (69.2) 


96 (80.7) 


9886 (87.4) 


Birth weight (gram) 












<2500 


156 (22.5) 


29 (8.2) 


34 (21.4) 


11 (9.2) 


471 (4.2) 


2500-3999 


477 (68.8) 


306 (86.9) 


1 1 9 (74.8) 


100 (84.0) 


1 0269 (90.8) 


=:4000 


60 (8.7) 


1 7 (4.8) 


6 (3.8) 


8 (6.7) 


567 (5.0) 


Delivery way 


Natural labor 


385 (55.6) 


223 (63.4) 


121 (76.1) 


70 (58.8) 


7327 (64.8) 


Cesarean 


295 (42.6) 


125 (35.56) 


36 (22.6) 


47 (39.5) 


3891 (34.4) 


Others 


13 (1.9) 


4(1.1) 


2 (1.3) 


2 (1.7) 


89 (0.8) 


Plurality 


Single 


668 (96.4) 


345 (98.0) 


158 (99.4) 


116 (97.5) 


11206 (99.1) 


Twins 


25 (3.6) 


7 (2.0) 


1 (0.6) 


3 (2.5) 


101 (0.9) 


Outcome 












Live birth 


625 (90.2) 


349 (99.1) 


133 (83.6) 


1 1 5 (96.6) 


1 1 272 (99.7) 


Still birth ' 


57 (8.2) 


3 (0.9) 


24 (15.1) 


3 (2.5) 


29 (0.3) 


Neonatal death 


11 (1.6) 


0 (0) 


2 (1.3) 


1 (0.8) 


6 (0) 



^CHD: congenital heart defect, ^ CL/P: cleft lip with or without palate, one case for unknown gender. 
*^Fetal deaths of at least 28 weeks' gestation, ^ Died within one week after birth. 
*P value <0.05, **P value <0.01 as compared with controls. 
doi:l 0.1 371 /journal.pone.0081 369.t001 



OR = 1.42, 95% CI 1.09-1.84). A two-fold increased risk of 
having a CL/P-affected offspring was observed in mothers older 
than 35 years old (crude OR = 2.14, 95% CI 1.30-3.52) compared 
with those of 25-29 years old, and adjustments for covariates 
slightly decreased the risk estimates (adjusted OR = 2.12, 95% CI 
1.26-3.57). There was no significant association between maternal 
age and equinovarus (for adjusted model: P=0.11). 

Maternal parity and risk of selected isolated birth defects 

As outlined in Table 3-4, the adjusted analysis revealed a 41 % 
and 58% decrease in the risk of delivering an infant with 
equinovarus among women giving birth to their second 
(OR =0.59, 95% CI 0.40-0.97), and third or more (OR = 0.42, 
95% CI 0.19-0.93) child respectively, as compared with primip- 
ara. Regarding CHD and Polydactyly, reduced risks were also 
observed in higher parity groups, whereas without statistical 
significances. The crude model for CL/P pointed to a 52% 
increased risk for women with parity &2 (OR= 1.52, 95% CI 
0.95-2.46), yet this risk estimate decreased and turned out to be 
non-statistical significant in the adjusted model, with an OR of 
1.34 (95% CI 0.81-2.23). 

Discussion 

In the current case-control study, we found that maternal age 
was significantly associated with CHD, Polydactyly and CL/P 



relevant pregnancy, and those with higher parity were less likely to 
have equinovarus-affected offsprings. 

The impacts of maternal age on birth defects have been long 
investigated, with inconsistent findings. An earlier study [16] 
conducted in the California population postulated that women 
39 years or more of age were twice as likely as 25-29 years old to 
have a child with isolated CL/P, to which our study results were 
similar. Later Reefhuis and Honein [17] had carried out a large 
sample-size study using a cohort of 1,050,616 singleton infants 
born from 1968 to 2000 in Atlanta, and found that young 
maternal age (14-19 years old) represented a significant risk factor 
for cleft lip (OR= 1.88, 95% CI 1.30-2.73). However, a recent 
published meta-analysis [18] on oral clefts extracted data from 13 
studied and suggested that, instead of early maternal age, mothers 
aged 40 years or over were 1.56 times more likely to have a 
newborn with non-syndromic CL/P (NSCL/P) compared with 
those aged between 20-29 years. 

Meanwhile, Reefhuis and Honein [17] also observed that 
advanced maternal age (35-40 years old) was associated with all 
heart defects (OR = 1.12, 95% CI 1.03-1.22). Similarly, Miller 
et al. [19] pointed an increased prevalence for several CHD 
phenotypes among women aging 35 years or older, with 
significant adjusted prevalence ratios ranged from 1.20 to 2.06. 
Yet our study did not replicate the above findings, on the contrary, 
we observed a reduced risk of CHD among younger mothers. 
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Table 2. Maternal factors of cases with Isolated birth defects and unaffected controls. 





Factors 


CHD ' 


Polydactyly 


CL/P " 


Equinovarus 


Control 




(n = 693, %) 


(n = 352, %) 


(n = 159, %) 


(n = 119, %) 


(n= 11307, %) 


Nationality 


Han 


634 (91.5) 


327 (92.9) 


147 (92.4) 


111 (93.3) 


10378 (91.8) 


Other 


59 (8.5) 


25 (7.1) 


12 (7.6) 


8 (6.7) 


927 (8.2) 


Household registry 


Local resident 


58 (8.4) 


13 (3.7) 


2 (1.3) 


5 (4.2) 


461 (4.1) 


Temporary population 


240 (34.6) 


100 (28.4) 


35 (22.0) 


25 (21.0) 


3867 (34.2) 


Floating population 


395 (57.0) 


239 (67.9) 


122 (76.7) 


89 (74.8) 


6979 (61.7) 


Antenatal examination (num) — 


0 


52 (7.5) 


27 (7.7) 


19 (11.9) 


4 (3.4) 


742 (6.6) 


1-5 


397 (57.3) 


221 (62.8) 


117 (73.6) 


75 (63.0) 


7010 (62.0) 


6-10 


213 (30.7) 


95 (27.0) 


19 (11.9) 


37 (31.1) 


3119 (27.6) 


>10 


31 (4.5) 


9 (2.6) 


4 (2.5) 


3 (2.5) 


436 (3.9) 


Pregnancy complication 


Gestational diabetes mellitus 


12 (1.7) ~ 


0(0) 


0(0) 


0(0) 


37 (0.3) 


Gestational hypertension 


13 (1.9) ' 


3 (0.9) 


0 (0) 


3 (2.5) 


96 (0.8) 


Mild preeclampsia 


4 (0.6) 


0(0) 


2 (1.3)* 


0 (0) 


23 (0.2) 


Severe preeclampsia 


12 (1.7) " 


8 (2.3) "" 


0(0) 


0 (0) 


68 (0.6) 



^CHD: congenital heart defect. 

^CL/P: cleft lip with or without palate. 

*P value <0.05, **P value <0.01 as compared with controls. 

doi:1 0.1 371 /journal.pone.0081 369.t002 



Unlike the widely explorations of CHD and CL/P, relatively litde 
was reported about the impacts of maternal age on Polydactyly and 
equinovarus. Our finding of a modest but significant association 
between young maternal age and a higher risk of Polydactyly was in 
accordance with one previous study [17] that expectant mothers 
under 20 years old was documented to at higher risk of carrying a 
Polydactyly oflspring (OR = 1.29, 95% CI 1.09-1.52). As for 
equinovarus, opposite conclusions had been drawn by two 



researches: Carey et al. [20] found that young maternal age 
(<20 years) was not associated with excess numbers of isolated 
talipes equinovarus cases while Parker et al. [21] observed 10% 
decreased risk in women younger than 23 years old in comparison 
with those of 23-34 years (Adjusted OR = 0.90, 95% CI 0.83-0.96). 
Nevertheless, our study did not replicate the above hypothesizes. 

Similarly, it had not been verified whether maternal parity could 
contribute to the occurrence of certain birth defects. Recently 



Table 3. Maternal age, parity and the risk of isolated CHD and Polydactyly. 



Factors CHD ' Polydactyly Control 





(n = 693, %) 


Crude OR 


Adjusted OR 


(n = 352, %) 


Crude OR 


Adjusted OR 


(n = 1 1 307, %) 






(95% CI) 


(95% CI) " 




(95% CI) 


(95% CI) " 




Age (year) _ ... 


<25 


159 (22.9) 


0.73 (0.60-0.88) 


0.73 (0.59-0.90) 


123 (34.9) 


1.43 (1.10-1.85) 


1.42 (1.09-1.84) 


4324 (38.2) 


25-29 


290 (41.8) 


1.00 (Reference) 


1.00 (Reference) 


114 (32.4) 


1 .00 (Reference) 


1 .00 (Reference) 


3550 (31.4) 


30-34 


164 (23.7) 


1.03 (0.84-1.25) 


1.03 (0.83-1.27) 


88 (25.0) 


1.40 (1.06-1.86) 


1.42 (1.07-1.88) 


2469 (21 .8) 


2:35 


80 (11.5) 


1.28 (0.99-1.66) 


1.29 (0.98-1.70) 


27 (7.7) 


1.10 (0.72-1.69) 


1.09 (0.71-1.66) 


964 (8.5) 


Parity (num) 


0 


390 (56.3) 


1.00 (Reference) 


1.00 (Reference) 


191 (54.3) 


1.00 (Reference) 


1.00 (Reference) 


5837 (51.6) 


1 


236 (34.1) 


0.80 (0.68-0.95) 


0.84 (0.71-1.01) 


130 (36.9) 


0.90 (0.72-1.13) 


0.88 (0.70-1.10) 


4403 (38.9) 


>2 


67 (9.7) 


0.94 (0.72-1.23) 


0.95 (0.71-1.28) 


31 (8.8) 


0.89 (0.60-1.31) 


0.81 (0.54-1.20) 


1 067 (9.4) 



^CHD: congenital heart defect. 

"^Adjusted for gestational age, birth weight, delivery way, plurality, outcome, household registry, gestational diabetes mellitus, gestational hypertension, severe 
preeclampsia, parity and age. 

■^Adjusted for gender, birth weight, severe preeclampsia and age. 
*P value <0.05, **P value <0.01 as compared with reference group. 
doi:l 0.1 371/journal.pone.0081 369.t003 
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Table 4. Maternal age, parity and the risk of isolated CL/P and Equinovarus. 





Factors 


CL/P ■ 






Equinovarus 






Control 


(n = 1 59, %) 


Crude OR 


Adjusted OR 


(n = 119, %) 


Crude OR 


Adjusted OR 


(n = 11307, %) 






(95% CI) 


(95% CI) " 




(95% CI) 


(95% CI) " 




Age (year) » 


<25 


61 (38.4) 


1.61 (1.12-2.35) 


1.29 (0.87-1.91) 


47 (30.3) 


1.73 (1.12-2.67) 


1.41 (0.90-2.22) 


4324 (38.2) 


25-29 


50 (31.4) 


l.CO (Reference) 


1.00 (Reference) 


36 (30.3) 


1.00 (Reference) 


1.00 (Reference) 


3550 (31.4) 


30-34 


25 (15.7) 


0.91 (0.56-1.47) 


0.89 (0.54-1 .46) 


22 (18.5) 


1.11 (0.65-1.89) 


1.25 (0.73-2.14) 


2469 (21 .8) 


2:35 


23 (14.5) 


2.14 (1.30-3.52) 


2.12 (1.26-3.57) 


14 (11.8) 


1.81 (0.97-3.37) 


2.13 (1.12-4.04) 


964 (8.5) 


Parity (num) 
















0 


79 (49.7) 


1.00 (Reference) 


1 .00 (Reference) 


75 (63.0) 


1 .00 (Reference) 


1.00 (Reference) 


5837 (51.6) 


1 


58 (36.5) 


0.97 (0.69-1.37) 


0.95 (0.66-1 .37) 


37 (31.1) 


0.65 (0.44-0.97) 


0.59 (0.40-0.89) 


4403 (38.9) 


>2 


22 (13.8) 


1 .52 (0.95-2.46) 


1.34 (0.81-2.23) 


7 (5.9) 


0.51 (0.24-1.11) 


0.42 (0.19-0.93) 


1 067 (9.4) 



^CL/P: cleft lip with or without palate. 

^Adjusted for gender, gestational age, birth weight, delivery way, outcome, household registry, antenatal examination, mild preeclampsia and age. 
*^Adjusted for birth weight, outcome, household registry, age and parity. 
*P value <0.05, **P value <0.01 as compared with reference group. 
doi:l 0.1 371 /journal.pone.0081 369.t004 



Duong et al. [22] pointed that nuUiparous mothers (those with no 
previous hve birth) were more likely to have infants with tetralogy of 
Fallot and sepals cardiac defects compared with mothers with one 
prior live birth, with significant ORs of 1.34 and 1.20. In contrast, 
Messer et al. [23] had proposed that slightly larger proportion of 
NSCL/P cases were born to women with three or more prior births 
in Texas during 1999-2003. However, our data did not support the 
relation between parity and the above mentioned defects, and only 
obtained a significant reduced risk for equinovarus in multipara, 
which was opposite to the viewpoint that an excessive number of 
parturitions may lead to reproductive organ damage and dysfunc- 
tion. Worth of note, the apparent protective effect of higher parity 
appeared to be of greater magnitude than that in previously report 
[21] which recruited 6,139 cases and 61,390 controls and showed 
significant adjusted ORs of 0.69 and 0.71 for birth order of second 
and third, respectively. 

We deemed that nonstandard criteria for the age and parity 
grouping might be partially responsible for the inconsistent findings. 
Generally, the risk estimates for maternal age were calculated using 
5-year categories with 25-29 years as the reference, while some 
studies introduced chflFerent classifications. For instance, Parker 
et al. [2 1] categorized the age into three groups with a wider range 
(<23, 23-34 and S:35 years old). The definitions of parity were also 
various. Currently, the term of parity refers to the numbers of 
maternal parturition processes of live birth [15]. However, the 
common argument for parity is probably because of the definition of 
live birth, more exactly speaking, the "surviving fetus" (or fetus with 
vigor) has no unified criteria under different biomedical and 
technological levels. In the US, for example, those whose gestational 
age >20 weeks or birth weight >500 gram were considered as 
"substantially viable fetus", while in the UK, parity is counted as the 
number of births beyond 24 weeks' gestation [24]. In China, this 
value was postponed to 28 weeks, which was also the lower limit 
value of perinatal period. As a consequence, there is, at the present 
time, no worldwide, uniform standard to define parity. In addition, 
selection bias from study arm, both cases and controls, could also be 
a possible explanation for discrepancies, since all the participants 
were recruited from different pools and therefore, population 
heterogeneities that were inherent in individual studies could not be 
ignored. 



Our study had several notable methodological advantages, 
namely: 1) the sample size of the present study was comparatively 
large and substantial to decrease the rate of false-positive or false- 
negative of the results, the 95% CIs for risk estimates were 
relatively narrow and suggested sound statistical precision; 2) for 
the widely coverage and objectively data collection from the 
population-based register, both cases and controls could be 
considered as representative of the general population in Baoan 
district. Moreover, the selection and information bias was 
rendered to minimum; 3) multivariable adjusted analyses were 
conducted to reduce the errors introduced by confounders, just as 
the liner regression analysis for the co-linearity problems. 

However, certain limitations were also needed to be acknowl- 
edged. Environmental conditions (such as occupational exposure of 
pollutants) or maternal certain personal behaviors (such as tobacco 
and alcohol consumption or nutritional deficiencies), as mentioned 
before, also had influence on the whole pregnancy outcome. 
Nevertheless, these information were not available for analysis, and 
it was difficult to derive such specific data because of the nature flaw 
of the Shenzhen Maternal and Child Health Management System, 
as the routine data collections were primarily used for administra- 
tive or governmental reports, most of the environmental and 
behaviors exposures were recorded as texts without uniform 
standard. Given the complexity of embryonic developments and 
the modest effect from single factor, it was possible that our findings 
were attributable to interactions among various components. Thus 
we stressed that the aforementioned results merited further 
comprehensive investigations. 

In summary, our findings indicated that maternal age and parity 
might play an important role in the development of certain isolated 
birth defects. Given the prevalence of these defects in Baoan district, 
our findings might help other researchers to conduct more 
etiological studies and plan more comprehensive strategies to 
reduce the occurrence of birth defects in further generation. 
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